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The IceCube Neutrino Observatory is a kilometer-scale detector currently under construction at the South 
Pole. In its final configuration the detector will comprise 5160 Digital Optical Modules (DOMs) deployed on 86 
strings between 1.5-2.5 km deep within the ice. While still incomplete, the detector has already recorded tens 
of billions of cosmic ray muons with a median energy of 20 TeV. This large sample has been used to study the 
arrival direction distribution of the cosmic rays. We report the observation of an anisotropy in the cosmic rays 
arrival direction at two different angular scales. The observed large scale anisotropy seems to be a continuation of 
similar structures observed in the Northern Sky by several experiments. IceCube observes also significant features 
on the angular scale of 20° - 30° that might be part of the larger scale structure. 



1. INTRODUCTION 

During the last twenty years several under- 
ground experiments have demonstrated the pres- 
ence, in the northern sky, of an anisotropy in the 
arrival direction of the cosmic rays (CRs) at en- 
ergies up to few hundred GeV p], [2]. Recently, 
both underground and surface array experiments 
have observed a sidereal anisotropic modulation 
on the order of 10 -4 in the multi-TeV range ( [3J, 
[1], [5] and [5J,). The origin of the anisotropy is 
unknown and observations of galactic cosmic ray 
anisotropy at different energy and angular scales 
have, therefore, the potential to reveal the con- 
nection between CRs and the postulated causes. 
Cosmic rays at these energies are almost entirely 
of Galactic origin and are expected to be isotropic 
due to the interaction with the Galactic magnetic 
field (GMF). Nevertheless, several theories have 
been postulated to explain the observations and 
the anisotropy can be imagined as a result of dif- 
ferent astrophysical phenomena. An anisotropy 
can be induced by the heliospheric magnetic field. 
At lower energies around a TeV, the heliosphere 
may be able to induce a CR excess in the direc- 



tion of the heliotail region (the so-called "tail-in" 
excess) [2]. At higher energies, the distribution 
of nearby recent supernova explosions has been 
postulated to be capable of creating a large-scale 
anisotropy [7] . In addition to these effects Comp- 
ton and Getting in 1935 8 predicted a dipole 
effect due to the motion of the Earth with re- 
spect to an isotropic CR plasma rest frame. Two 
anisotropies can be postulated, one due to the 
Earth's motion around the Sun (Solar Dipole ef- 
fect), that would appear in solar time, and the 
other coming from the motion of the solar sys- 
tem around the galactic center (Galactic Comp- 
ton Getting effect). This second effect would ap- 
pear in sidereal time (ST) and it would be maxi- 
mum with the CRs being at rest with respect to 
the galactic center. 

In this contribution we present the observations 
of the CR anisotropy by IceCube at two differ- 
ent angular scales. With this measurement we 
produced the first skymap in cosmic rays of the 
southern sky. 
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2. DATASET 

The IceCube Neutrino Observatory currently 
under construction at the South Pole will be com- 
pleted at the beginning of 2011 with 86 strings 
and 5,160 DOMs. While IceCube is designed to 
observe the sources of extragalactic neutrinos, it 
is also sensitive to TeV downward-going muons 
produced by cosmic ray air showers. Since the 
start of data taking in 2007, IceCube has recorded 
tens of billions of muon events, and has accu- 
mulated the largest sample of TeV cosmic rays 
ever recorded in the southern hemisphere. The 
results presented here have been obtained analyz- 
ing the data collected in two periods, from June 
2007 to March 2008 and from May 2008 to April 
2009, in which IceCube was operating with 22 
and 40 strings, respectively. The events used in 
the analysis have been reconstructed by an online 
likelihood-based reconstruction algorithm. The 
event rate is ~ 240 Hz for IceCube-22 and ~ 750 
Hz for IccCube-40, with a median cosmic ray en- 
ergy of about 20 TeV and a median angular res- 
olution of 3° , where this value represent IceCube 
angular resolution for atmospheric muons with no 
quality cuts applied and has not to be confused 
with the angular resolution for neutrinos (better 
than 1° for IceCube-40). 

These large statistics data sets allow for sensitive 
studies of the arrival direction distribution of the 
cosmic rays. 

3. DATA ANALYSIS AND RESULTS 

3.1. Large scale anisotropy 

IceCube recently reported the first observation 
of the large angular scale anisotropy in arrival di- 
rection of CRs for the southern hemisphere [9]. 
The analysis determines the relative intensity of 
the cosmic ray arrival direction in each declina- 
tion band and it has to account for the spurious 
effects derived from the detector exposure asym- 
metries and non-uniform time coverage as well as 
for diurnal and seasonal variations of atmospheric 
conditions. Fortunately, due to its favoured po- 
sition at the South Pole, IceCube measurements 
of the sidereal variation is not affected by diurnal 
modulations because the whole sky is fully visi- 



ble to the detector at any given time and because 
there is only one day and one night per year; on 
the other hand, the seasonal variation in the cos- 
mic ray event rate is slow and does not affect the 
daily muon intensity significantly. 
The remaining spurious effects which must be ac- 
counted for are an asymmetry in the detector ge- 
ometry and non-uniformity in the time coverage 
of the data. The combinations of these two ef- 
fects might induce an azimuthal asymmetry and 
mimic a sidereal anisotropy. This effect is cor- 
rected normalizing the azimuthal distribution by 
re-weighting each event with the ratio between 
the average number of events and the number 
of events in the corresponding local azimuthal 
bin. Moreover, since the local azimuth distribu- 
tion varies with the zenith angle, the sky is di- 
vided into four zenith bands with approximately 
the same number of events per band. The weight- 
ing is applied within each band to remove the de- 
tector asymmetry. 

Figure [T] shows the relative intensity in arrival 
direction of the cosmic rays: on the top the rela- 
tive intensity map obtained from the ~ 4.6 x 10 9 
events collected by IceCube-22 [9] and on the 
bottom the preliminary map obtained from the 
~ 12 x 10 9 events collected by IceCube-40. Since 
the arrival direction distribution is dominated by 
the zenith angle dependence of the flux, the maps 
have been obtained by normalizing to unity each 
declination belt of 3°, which corresponds to the 
angular resolution of the data. The two maps 
show the same anisotropy features and they both 
appear to be a continuation of the observed mod- 
ulation in the northern hemisphere. To quan- 
tify the scale of the anisotropy the right ascen- 
sion distribution has been fitted with a first and 
second-order harmonic function in the form of 
Egfj^Aj x cos(i((a) - 0,-))) +B, where A, is the 
amplitude, <f>i is the phase and B is a constant. 
The fit parameters as well as the quality of the fit 
for IceCube-22 are reported in [5]. 

Figure [2] shows the right ascension modulation 
of the relative intensity in arrival direction of cos- 
mic rays for IceCube-22 (black dots) [9] as well 
as the preliminary result for IceCube-40. The er- 
ror bars represent the statistical uncertainties and 
the lines are the fit to the data. The gray band in- 
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Figure 1. Sky-map of the relative intensity in ar- 
rival direction of cosmic rays for IceCube-22 (top) 
[9] , and preliminary sky- map of the relative inten- 
sity for IceCube-40 (bottom), in equatorial coor- 
dinates. A gaussian smoothing has been applied 
to the map for visualization purposes. Note that 
since the declination belts in the equatorial map 
are treated independently, the maps provide only 
information on the relative modulation of the ar- 
rival direction of cosmic rays along the right as- 
cension. 



dicates the systematic uncertainties for IceCube- 
22 estimated taking into account daily and sea- 
sonal variations and non-uniform time coverage 
in the data. 

3.2. The Solar Dipole effect 

Any observer moving through a plasma of 
isotropic cosmic rays should observe a distinct 
difference in intensity between the direction of 
the velocity vector and the opposite direction. In 
1935 Compton and Getting predicted that the in- 
tensity of CRs should be observed like a dipole 
anisotropy [8] ; in this model an excess flux should 
appear with a maximum in right ascension be- 
tween 290° and 340° and a minimum in right as- 
cension between 110° and 160° [3]. The large 
scale anisotropy observed in IceCube (FigEJ) can 
not be described as a pure dipole and the ob- 
served excess is not in the direction of motion of 
the solar system around the galaxy. We can con- 
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Figure 2. Right ascension modulation (for 
-72° < 5 < -38°) of the relative intensity 
in arrival direction of cosmic rays for IceCube- 
22 (black dots) 9 and preliminary result for 
IceCube-40 (red dots). The data are shown with 
statistical uncertainties, and the lines represent 
the fit to the data. The gray band indicates the 
systematic uncertainties for IceCube-22. 



elude that the galactic Compton-Gctting effect is 
at most one of several contributions to the side- 
real anisotropy. 

We also expect a dipole anisotropy caused by the 
Earth's motion around the Sun. The expected 
anisotropy is of order 10~ 4 . If we represent the 
relative intensity as a function of the angular dis- 
tance from the Sun in right ascension (defined as 
the difference between the right ascension of each 
event and the Sun position), we expect, over a full 
year, an excess in the direction of motion of the 
Earth around the Sun (at 270°) and a minimum 
in the opposite direction. Figure [3] shows the 
modulation observed in IceCube-40. The black 
dots represent the experimental data with the 
corresponding statistical uncertainties. The data 
have been fitted with a first-order harmonic func- 
tion and the result of the fit is represented by the 
black line. The red line shows the expected solar 
dipole effect. In the plot the Sun is located at 
0°and the velocity vector of the Earth is at 270° 
which is where the maximum of the modulation 
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Figure 3. Preliminary modulation of the cosmic 
ray relative intensity as a function of the rela- 
tive arrival angle from the Sun in right ascension 
observed in IceCube-40. The black dots are the 
experimental data with the corresponding statis- 
tical error bars and the line represent the fit to 
the data. The red line represents the expected 
modulation due to Earth's motion. 



is observed, as expected. 

The dipole effect caused by Earth's motion 
around the Sun is visible only when the arrival 
directions are plotted in a frame where the Sun 
position is fixed in the sky. A fixed signal in this 
coordinate system is completely washed out in 
sidereal time over the course of one year. There- 
fore, the solar dipole effect is not visible in the 
sidereal skymap. 

3.3. Medium scale anisotropy 

Recent results from gamma ray experiments 
have shown evidence of intermediate (~ 10°— 30°) 
scale structure in the CR arrival direction distri- 
bution. Milagro observed two localized regions 
with flux excess of significance > 10a in the to- 
tal data set of 2.2 x 10 11 events recorded over 
7 years. The "hot" regions have fractional ex- 
cesses of the order of several times 10 -4 relative 
to the background [10 . The same structures have 
been observed by ARGO-YBJ 11 and Tibet III 



[j"2] . The origin of these excesses is still unknown 
but there seems to be some indications that the 
medium scale features might be part of the larger 
scale structure. 

With the same technique used in the gamma ray 
field, IceCube can search for similar excess in the 
southern sky. The analysis technique is based on 
the background estimation from real data. The 
method can have different implementations and 
its general features are described in [T3]. The 
analysis used in IceCube can be summarized in 
few steps: 

• a signal map is created from the arrival di- 
rections of the events; 

• a background map is calculated from the 
real data by randomly assigning detected 
event times to local arrival directions 
(time scrambling). The uncertainty in the 
background estimation is reduced by re- 
sampling the event time 20 times. 

• the statistical significance of the result is 
calculated using the method of Li and Ma 

The length of the time window used to scramble 
the time of the events determines the maximum 
angular scale we can probe. Since we are inter- 
ested in searching for features of angular scale 
~ 10° — 30° we choose a time window of 4 hours. 
In this way the time window is used to filter the 
signal coming from the large scale structure since 
only the anisotropy features with angular scale 
smaller than 60° are visible. Moreover, in order to 
maximize the sensitivity, signal and background 
maps are smoothed on scales that correspond to 
the feature size. 

Figure |4] shows the preliminary significance 
skymap in equatorial coordinates for a time win- 
dow of 4 hours obtained with ~ 1.9 x 10 10 events 
collected in one year of IceCube-40. The map 
has been smoothed using an optimal radius of 25° 
(determined by examining the data) to improve 
the sensitivity to large features. The smoothing 
is performed adding counts from all pixels within 
a certain radius; due to this procedure the pixels 
in the map are strongly correlated. The skymap 
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Figure 4. Preliminary significance skymap in 
equatorial coordinates for a time window of 4 
hours obtained with IceCube-40 data. A 25° 
smoothing has been applied to improve the sen- 
sitivity to large features. The skymap shows a 
broad excess region around right ascension 120°, 
with an equally strong deficit around right ascen- 
sion 220°. 



shows a broad excess region around right ascen- 
sion 120°, with an equally strong deficit around 
right ascension 220°. 

Figure [5] shows the preliminary RA projection 
of the number of events for the signal and the 
expected background for —45° < S < —30°. The 
plot is made using independent 15°<5 x 1° RA bins 
(i.e. no smoothing applied). The data used for 
this figure has been chosen to include only full 
days in order to achieve an approximately uni- 
form exposure as a function of RA. As can be 
seen, the background estimate as calculated via 
the time scrambling technique agrees well with 
the data. The excess and the deficit region shown 
in the skymap are visible in the raw data. 

4. CONCLUSIONS 

IceCube observes, for the first time in the 
southern hemisphere, an anisotropy in the arrival 
direction of the galactic cosmic rays. The large 
scale structure seems to be a continuation of a 
previously observed modulation in the northern 
hemisphere. The IceCube skymap also reveals the 
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Figure 5. Preliminary right ascension projection 
of signal and background events for —45° < 6 < 
—30°. The plot is made using independent 15° 6 x 
1° RA bins (i.e. no smoothing applied). The 
data used here include only the full days to ensure 
a uniform exposure in RA. The two significant 
regions are already visible in the raw data. 



presence of an anisotropy in intermediate scale 
(~ 10° — 30°) with a broad excess region around 
right ascension 120°, and an equally strong deficit 
around right ascension 220°. The origin of these 
anisotropics is unknown. However there might be 
multiple superimposed causes, depending on the 
cosmic ray energy and the angular scale of the 
anisotropy. The postulated Compton-Getting ef- 
fect [5] should give rise to an energy independent 
dipole anisotropy with the maximum in the direc- 
tion of the solar system and the deficit in the op- 
posite direction. The measured modulation can 
not be described as a pure dipole and the observed 
excess is not in the direction of motion of the solar 
system around the galaxy; therefore, we can con- 
clude that the galactic Compton-Getting effect, 
if present at all, is overshadowed by other effects. 
One could also imagine a scenario in which the 
cosmic ray excess might be associated with dif- 
fuse particle flows coming from nearby cosmic ray 
sources like Vela SNR. Recently the MILAGRO 
excess has been interpreted by invoking Geminga 
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pulsar as a possible source [IS], [IS]- However, 
it is likely that a localized excess of multi-TeV 
CRs is originated at close distances. In Lazarian 
& Desiati it is proposed that both the localized 
excess observed by Milagro toward the direction 
of the heliotail (i.e. the portion of the heliosphere 
downstream the interstellar wind) and the broad 
excess of cosmic rays observed at sub-TeV ener- 
gies from the same direction (tail-in excess) [2] 
are generated by first-order Fermi acceleration via 
stochastic magnetic reconnection in the heliotail 
up to about 10 TeV [T7]. Reconnection is gener- 
ated when solar magnetic polarity reversal regions 
due to the 11-year solar cycles are compressed by 
the solar wind in the heliotail. In the reconnec- 
tion regions, magnetic mirrors form where accel- 
eration might be efficient to distort spectrum and 
arrival direction of the cosmic rays. 
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